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Controlled doping of SrCu2(B03)2, a faithful realization of the Heisenberg spin-1/2 antiferromag- 
net on the Shastry-Sutherland lattice, with non-magnetic impurities generates bound-states below 
the spin gap. These bound-states and their symmetry properties are investigated by exact diago- 
nalisation of small clusters and within a simple effective model describing a spinon submitted to an 
attractive extended potential. It is shown that Raman spectroscopy is a unique technique to probe 
these bound-states. Quantitative theoretical Raman spectra are numerically obtained. 

PACS numbers: 75.10.Jm, 75.40.Gb 
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I. INTRODUCTION AND MOTIVATIONS 

SrCu2 (603)2 is an experimental realization of a spin 
1/2 antiferromagnet living on the two-dimensional (2d) 
Shastry-Sutherland lattice (SSL). Its properties can be 
well described with the Heisenberg model: 
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where J (resp. J') is the exchange within (resp. be- 
tween) dimers. The experimental compound has a ratio 
a = J' /J slightly above^ 0.6. The SSL has been first 
introduced theoretically as an example of a 2d antiferro- 
magnet whose ground-state for small enough a is exactly 
knowi*i: It is simply the product of singlets on each dimer 
up to^ a ~ 0.7. Indeed, experiments on SrCu2 (603)2 in- 
dicate a finite spin gap. For larger a, a more usual 2d 
Neel-ordcred phase is stabilized, possibly with an inter- 
mediate plaquette phase. Previous Raman experiments 
on this compound^ indicate that some structure in the 
spectrum can only be explained by taking into account 
a realistic model, including for instance Dzyaloshinski- 
Moriya (DM) interactions^. However, these DM terms 
are too small to have any sizeable effects on the low- 
energy Raman spectra^ discussed here and thus will be 
neglected hereafter. Our goal will be to show that, within 
a simple Heisenberg model, doping non-magnetic impu- 
rities generates new spectroscopic signatures below the 
2-magnon continuum whose features should persist in a 
more realistic model.— 

SrCu2 (603)2 is a Mott insulator and its doping with 
mobile carriers is predicted to lead to superconductivity^. 
In this paper, we shall rather consider the controlled dop- 
ing of SrCu2 (603)2 with non-magnetic static impurities 
such as Zinc or Magnesium atoms substituted for Copper 
atoms on a very small fraction of the N lattice sites. Such 
atoms acting as vacant sites provide accurate local real- 
space probes (which can be considered as independent 
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FIG. 1: (Color online) Shastry-Sutherland lattice. ± signs 
correspond to Raman coupling in the {a'b') polarization case. 



for low enough impurity concentration). As known in 
strongly correlated systems, a small doping of the parent 
compound can bring crucial informations about its in- 
trinsic properties^. Substituting Cu with non-magnetic 
Mg impurities has been recently performe d^^'^"^ . Neutron 
scattering experimentsii have then shown the appear- 
ance of new magnetic excitations into the singlet-triplet 
gap. Here, we shall focus on Raman spectroscopic tech- 
niques which, by probing AS" = excitations using light 
scattering, offer a unique way to identify the local re- 
sponse to a doped impurity. In particular, we will show 
that transitions between different S=l/2 bound-states 
can be identified in the Raman spectra. Such features can 
be interpreted within a simple phenomenological model 
describing the attractive potential between a liberated 
spinon and the impurity (vacant site). 



II. THEORETICAL FRAMEWORK 

A. Raman scattering 

The theory of Raman light scattering is most simple 
when the photon energy is much smaller than the Mott 
gap. In such a case, it is legitimate to use the Loudon- 
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Fleury approximation and the Raman operator reads: 
7?. = ^ 7(em • dy )(eo„t • dy ) Si ■ Sj 

nn 

+ X! 'y'i^^n ■ d.y)(eo«t ■ dy ) S,; ■ Sj (2) 



where ei„ and Gout are the polarization vectors of the 
incoming and scattered Hght, and d^j is the unit vector 
connecting two sites i and j. This operator only cou- 
ples to zero-momentum singlet excitations. In principle, 
the coupling constants 7 and 7' could depend on the ex- 
change values^, although such a dependance is neglected 
in most studies. However, some simplifications occur for 
an {a'b') polarization (see Fig. [T|): indeed, in this case, 
the geometry of the compound implies that only the J' 
bonds contribute to the Raman operator. Moreover, the 
dominant intensity occurs when triplets can be created on 
dimer bonds which, in contrast, is not allowed in (ab) po- 
larization^ii. In the following, we will restrict to this po- 
larization where all Raman coupling constants are equal 
up to a sign (given in Fig. [T]). 

At zero temperature, the Raman intensity is given by 
the dynamical correlations of the Raman operator: 
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where j^E'o) is the ground-state (GS) of the system and 
the sum runs over the excited states \n) with energy E„. 
Note that because of the symmetry of the Raman opera- 
tor, only singlet states with zero momentum contribute. 
Moreover, in the chosen polarization, only states which 
are odd with respect to reflections along a' or 6' axis give 
a signal. 
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B. Results for the pure compound 

As studied in Ref. jj, the Raman spectrum of the un- 
doped material shows 4 sharp peaks at 1.25, 1.9, 2.3 and 
2.9 times the spin gap Agi. Naively, one would expect 
the Raman spectrum to start at the 2-magnon continuum 
(i.e. twice the spin gap) since it is a singlet operator. In 
fact, a singlet bound state made of 2 triplets does ex- 
ist on this lattice^ and was identified as the low-energy 
stated. The first two peaks are attributed to two-triplet 
bound states, while the 2 others are interpreted in terms 
of three-particle excitations. Although the high-intensity 
peaks are above the 2-magnon continuum (starting at 
2A01), the very small magnon dispersion leads to a very 
narrow continuum so that the two higher-energy peaks 
are observable. 

In Fig. [21 we present our Exact Diagonalizations spec- 
tra for various a. As expected, for small a, the finite-size 
effects are rather weak due to the large energy scales (i.e. 



FIG. 2: (Color online) Raman spectra for pure SSL for various 
Of on TV = 16 and 32 clusters. The spin gap Aoi (~ 0.95J, 
0.80 J, 0.50 J in (a), (b) and (c), respectively) is used as a unit 
of frequency. An artificial width e = 0.01 has been given to 
the delta peaks. 



short correlation lengths), so that results are almost iden- 
tical on = 16 and A^ = 32 clusters. Clearly, spectral 
weight is present below the 2-magnon threshold. On the 
other hand, for larger a, finite-size effects become size- 
able so that a direct comparison with experimental values 
is difficult. However, one can still notice spectral weight 
well below the 2-magnon continuum, corresponding to 
the above mentionned singlet bound states. Note also 
that, for A^ = 32, the first peak is located at an energy 
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1.2 Aqi as seen in experiments. 



C. Results for the doped case 

We now turn to the doped case for which additional 
low-energy states appear. We shall assume a small 
enough impurity concentration so that a single impu- 
rity description becomes legitimate. Using a varia- 
tional approach, El Shawish and BoncaM have proposed 
anisotropic spin-polaronic states with a finite spatial ex- 
tension around the impurity. Adding a single impurity 
(i.e. creating a vacant site by removing a spin) indeed 
generates a polarization 1/2 that will distribute around 
the impurity. On Fig. [3l we show the GS magnetiza- 
tion pattern computed exactly on a = 32 cluster with 
one impurity. The local polarization oscillates from one 
site to the other similar to the Friedel oscillations re- 
ported around a localized triplet in the 1/8 plateaui^. 
Note however that each strong dimer has a global pos- 
itive polarization. Our results are in good agreemen^i^ 
with the GS variational estimate of Ref. 0- These au- 
thors have also found low-energy states appearing below 
the undoped spin gap. Since these states are well local- 
ized, we expect that finite-size calculations can provide 
accurate results for these excitations as well. In particu- 
lar, we expect a much better accuracy than in the pure 
case. 
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FIG. 3; (Color online) Magnetization pattern around the im- 
purity (marked as x) on iV = 32 SSL for a = 0.6 (Sf = 1/2). 
Periodic boundary conditions are used. Since the reflection 
around the b' axis is still a good symmetry, we observe iden- 
tical values on both sides. 



sensitive to AS* = 1 transitions); moreover, since these 
states can be connected by AS* = transitions, they are 
also expected to be Raman active, which is the main 
purpose of our study. Fig. [¥] shows the Raman spectra 
obtained by Exact Diagonalizations on A = 16 and 32 
clusters with one impurity for various a. The spectra are 
plotted as a function of w/Aqi, where the spin gap of the 
pure sample Aqi is very close to the spin gap Ai 3 in the 
doped cluster (both computed on the largest cluster). In 
the small a regime, all energy scales are well separated 
and finite-size effects are negligible so that one can un- 
derstand all features. For the pure system, as discussed 
above, there is a bound-state below 2Aoi and then an- 
other peak at 2Aoi corresponding to a singlet excitation 
made of two distant triplets. Since the triplets have a 
small dispersion, there is no continuum above. In the 
presence of a single impurity, these features still repre- 
sent much of the spectral weight but, in addition, new 
peaks appear below the spin gap. Namely, one can make 
a singlet excitation by creating a triplet in the bulk while 
flipping the polarization cloud (so that the number of 
such peaks scales as the number of dimers A^ = ^ !)■ 
Moreover, there is a possibility to form a bound-states of 
these two excitations, which can even lower the energy 
below the spin gap (see Fig.[4|), as found numerically. An- 
other feature of this additional spectral weight is that it 
scales with the concentration of impurity, i.e. is reduced 
by 2 when doubling the system size. 

In order to be more quantitative, we compare these 
bound-state energies with the variational results of 
Ref. dil for a 0.62 where the two lowest 5=1/2 



states have excitation energies 0.238J and 0.264J, while 
the undoped spin gap is 0.450 J. In our exact calculations 
on A^ = 32 cluster for the same a, we find that the three 
lowest excitations have a total spin 1/2 and are located 
at 0.217, 0.245 and 0.268 J, while the spin gap is 0.479J 
in good agreement with these variational results, except 
that we have an additional low-energy state. Concerning 
their Raman signatures, one has to discuss their symme- 
try properties: in the presence of a single impurity, the 
translation symmetry is lost and only one reflection along 
dimers (i.e. along a' or b') remain. Since the Raman op- 
erator is odd with respect to this reflection, a simple in- 
spection at the odd/even character of these excited states 
allow to determine if they are/are not Raman active. In 
Fig.mja-b); we have indicated all even and odd Stot = 1/2 
states below the first 5*40* = 3/2 state on TV = 16 and 32 
clusters. A careful comparison with the Raman spectra 
does confirm that only states with Stot = 1/2 and which 
are odd w.r.t. reflection give Raman peaks at low energy. 



III. EFFECTIVE MODEL FOR THE 
BOUND-STATES 



The existence of several Sz = 1/2 low-energy states has 
led to various experimental signatures: by flipping spins, 
they can be observed in neutron experiment^i (which is 



Interestingly, a phcnomenological description of these 
bound-states can be given similarly to the case of doped 
quasi-lD CuGeOgi^iil. For J' = 0, a free 5 = 1/2 
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FIG. 4; (Color online) Raman spectra for SSL doped with a 
single non-magnetic impurity for various q on = 16 and 
32 clusters. The spin gap Aoi (~ 0.95 J, 0.80 J, 0.50 J in (a), 
(b) and (c), respectively) is used as a unit of frequency. An 
artificial width e — 0.01 has been given to the delta peaks. 



(spinon) is located next to the vacant site, on the broken 
dimer bond. Switching on J' allows this spinon to delo- 
calize with a hopping term of order (J')^/4J. However, 
each time the spinon hops one dimer away from the im- 
purity, a strong bond is broken resulting in an additional 
"string" energy cost (J — J'). Therefore, the physics is 
similar to a particle in a linear potential and bound-states 
can occur. Of course, when the string energy exceeds the 



spin gap. the whole picture breaks down and the spinon 
can "escape" in a flat potential by the spontaneous cre- 
ation of a spinon-antispinon pair out of the vacuum. We 
have considered this effective quantum mechanical model 
on an effective square lattice for one particle allowed to 
hop with amplitude a^/A on its neighboring sites (ex- 
cept between the impurity dimer and the neighboring 
dimer facing the vacant site), and with a potential en- 
ergy equals to V{r) = J min((l — a)d(r), 1) where d is 
the Manhattan distance from the impurity dimer. This 
one-particle problem can be easily solved on large clusters 
and its spectrum is presented on Fig.[S{c). For vanishing 
J', the low-energy S=l/2 spectrum is extensively degen- 
erate: the Nd = Y ~ ^ triplets of excitation energy J 
located on the remaining bonds can be combined with 
the impurity spin in total S=l/2 states (degenerate with 
their S=3/2 counter-parts not described by the model). 
The next set of states which appears at energy 2J (cor- 
responding to 2 isolated triplets) and above are also not 
described by the effective model. Switching on J' lifts 
the degeneracy of the first group of Nd S=l/2 states re- 
sulting in a rich spectrum well described by the effective 
model for which the N —>■ oo limit can be taken. Com- 
mon features are observed both in the microscopic and 
the effective models such as (a) bound-states due to the 
short-range string-like part of the potential and (b) a con- 
tinuum of S=l/2 excitations above the spin gap (given 
in Fig. [Sja-b) by the lowest S=3/2 state). In fact, the es- 
timation of the spin gap from high-order perturbation^^, 
Aoi ~ J (l — — — 1'^"') ' agi'ccs very well with the 
ED value obtained on the largest = 32 cluster up to 
a = 0.6 and with the effective model up to a = 0.3. For 
small a, the number of bound states below the spin gap is 
four: in the effective language, it corresponds to a spinon 
delocalized on one of the nearest-neighbor dimers of the 
impurity dimer. The situation for larger a is less clear as 
many states go down, possibly with stronger finite-size 
effects on ED data. Still, our numerical data would be 
compatible with up to 12 bound-states, some of them be- 
ing very close to the spin gap. However, with our choice 
of polarizations and due to the selection rule, only odd 
states are Raman active, which gives 2 (resp. 6) low- 
energy states for a 0.2 (resp. a ^ 0.6). 



IV. FINITE IMPURITY CONCENTRATION 

To finish, we quickly address the case of a finite im- 
purity concentration for which a two-impurity effective 
interaction starts to operate. Since each impurity will 
create a spin 1/2 polarization around it, it is natural 
to expect some ordering, possibly at low temperature. 
For instance, antiferromagnetic (AF) ordering occurs in 
doped CuGeOji^ and has been predicted for CaV409^. 
One can imagine an effective diluted spin 1/2 model with 
a very small exchange interaction (typically, the overlap 
between two polarization clouds is exponentially small 
when the impurities are quite far), and such a model is 
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FIG. 5: (Color online) Low-energy S=l/2 excitations vs a. 
The dashed line indicates the perturbative estimation of the 
spin gap (see text), (a-b) Microscopic model on = 16 and 
32 SSL. The states are classified according to their refiection 
symmetry. The two lowest S=3/2 states also shown give the 
onset of the continuum, (c) Effective model (see inset for a 
schematic plot of the effective potential) on a TV = 625 square 
lattice of dimers for which finite-size effects are not visible. 



fcctivc exchange can be obtained by considering various 
configurations of two impurities. Our results (data not 
shown) indicate that, when both impurities are located 
on vertical dimers, the effective interaction decreases very 
fast and is mostly AF. Since these bonds belong to the 
same sublatticc, the presence of strong frustration can 
prevent magnetic ordering. On the contrary, when both 
impurities are located on different dimcr types (one verti- 
cal and one horizontal) , the effective interaction is mostly 
ferromagnetic, thus competing with the AF and often re- 
sulting in a disordered state. Anyhow, if magnetic order 
occurs, it would be at a temperature much below the ef- 
fective energy scale, which is already quite small. This 
argument can be generalized to an arbitrary distribution 
of impurities. In conclusion, we predict the absence of 
magnetic ordering, even for quite large doping and at ex- 
tremely low temperature, which seems compatible with 
experiments^^. 



expected to order at low-temperature. 

In order to estimate the effective coupling constant 
between two impurities at site r and r', we start from 
the formula Jeff = 2{E{]]) - E{]i)). Now, E{]]) = 
'^{ij) Jiji^i ' "^i)' where the expectation value is calcu- 
lated in the triplet state with two impurities. Choosing 
the polarization along z, and in order to get an estimate, 
we can approximate the expectation value by {Sf){S^). 
We further approximate the polarization of a given site as 
the sum of the polarizations coming from the two impu- 
rities, which would be true to first order in perturbation, 
leading to 

{Sls+St,^{Sl,+S^^^ = (4) 
{Sl,){Sl,) + {Sta){Sl^) + {Sl,){Sl^) + {Sl^){Sl,)^ 

where (S'f i) {{8^2)) ^he average magnetization cre- 
ated at site i by a single impurity located at r (r'). 
The same appromixation for the correlation functions of 
i?(||) leads to the same terms with the same sign for 
the 1-1 and 2-2 terms, and opposite signs for the 1-2 and 
2-1 terms. The 1-1 and 2-2 terms will drop from the 
difference, leading to the formula: 

Jeff(r, r') = 2Y,M{Sl,){Sl2) + (5- 2)(^|,i)) (5) 
From the polarization shown in Fig. [31 the resulting ef- 



V. CONCLUSION 

In summary, doping a Shastry-Sutherland lattice with 
non-magnetic impurities leads to novel low-energy states 
below the spin gap, that could be probed by Raman spec- 
troscopy. We have proposed a simple effective model to 
understand bound-state formation as binding of a spinon 
to an impurity site. For a particular polarization of light, 
some of these bound-states are Raman active with size- 
able spectral weights: 2 of these bound-states are located 
well below the continuum (and exist for any J'/ J), while 
up to 4 more bound-states could be observed for realis- 
tic parameters. Moreover, since the effective interactions 
between impurities is frustrated, we expect no magnetic 
ordering down to extremely low temperature. 



Acknowledgements 

We acknowledge B. Gaulin for very useful discussions. 
We thank CALMIP (Toulouse) and IDRIS (Paris) for 
allocation of cpu time. This work has been supported by 
the French Research Agency (ANR), the Swiss National 
Fund, by MaNEP, and the region Midi-Pyrenees through 
its Chaire d 'excellence Pierre de Fermat. FM also thanks 
LPT (Toulouse) for hospitality. 



* Electronic address: |capponi@irsamc.ups-tlse.fr] 
^ B. S. Shastry and B. Sutherland, Physica B & C 108, 1069 
(1981). 

^ S. Miyahara and K. Ueda, J. Phys.: Condens. Matter 15, 
R327 (2003) and references therein. 



^ A. Gozar, B. S. Dennis, H. Kageyama, and G. Blumberg, 
Phys. Rev. B 72, 064405 (2005). 

* Y. F. Cheng, O. Cepas, P. W. Leung, and T. Ziman, Phys. 
Rev. B 75, 144422 (2007); V. V. Mazurenko, S. L. Sko- 
rnyakov, V. I. Anisimov, and F. Mila, Phys. Rev. B 78, 



6 



195110 (2008). 

^ In the presence of DM interactions, S > spin states might 
become Raman active. However, we are not concerned here 
with such excitations which (i) lie at higher energies than 
the energy range considered here and (ii) would only carry 
a small spectral weight. In addition, (iii) a perturbative 
DM interaction has an even weaker effect on the Raman- 
active low-energy singlet-like excitations. 
P. Lemmens, M. Grove, M. Fischer, G. Giintherodt, V. N. 
Kotov, H. Kageyama, K. Onizuka, and Y. Ueda, Phys. 
Rev. Lett. 85, 2605 (2000). 

C. Knetter, A. Biihler, E. Miiller-Hartmann, and G. Uhrig, 
Phys. Rev. Lett. 85, 3958 (2000). 
* P. J. Freitas and R. R. P. Singh, Phys. Rev. B 62, 14113 
(2000). 

® B.S. Shastry and B. Kumar, Prog. Theor. Phys. Suppl. 
145, 1 (2002); J. Liu, N. Trivedi, Y. Lee, B. N. Harmon, 
and J. Schmalian, Phys. Rev. Lett. 99, 227003 (2007). 
S. Sachdev and M. Vojta, Phys. Rev. B 68, 064419 (2003). 
S. Haravifard, S. R. Dunsiger, S. El Shawish, B. D. Gaulin, 
H. A. Dabkowska, M. T. F. Telling, T. G. Perring, and J. 
Bonca, Phys. Rev. Lett. 97, 247206 (2006). 
K. Kodama, M. Takigawa, M. Horvatic, C. Berthier, H. 
Kageyama, Y. Ueda, S. Miyahara, F. Becca, and F. Mila, 



Science 298, 395 (2002). 

A. A. Aczel, G. J. MacDougall, J. A. Rodriguez, G. M. 
Luke, P. L. Russo, A. T. Savici, Y. J. Uemura, H. A. 
Dabkowska, C. R. Wiebe, J. A. Janik, and H. Kageyama, 
Phys. Rev. B 76, 214427 (2007). 

S. El Shawish and J. Bonca, Phys. Rev. B 74, 174420 
(2006). 

For a = 0.62, the relative error between ground-state en- 
ergies with a single impurity obtained with ED on A'^ = 32 
or in the variational approach of Ref. [3 is 1.4%. 
" G. Els, G. S. Uhrig, P. Lemmens, H. Vonberg, P. H. M. 
van Loosdrecht, G. Guentherodt, O. Fujita, J. Akimitsu, 
G. Dhalenne, A. Revcolevschi, Europhys. Lett. 43, 463 

(1998) . 

^""^ D. Augier, E. S0rensen, J. Riera, and D. Poilblanc, Phys. 
Rev. B 60, 1075 (1999). 

S. Miyahara and K. Ueda, Phys. Rev. Lett. 82, 3701 

(1999) . 

" J.-P. Renard, K. Le Dang, P. Veillet, G. Dhalenne, A. 
Revcolevschi, and L.-P. Regnault, Europhys. Lett. 30, 475 
(1995). 

^° S. Wessel, B. Normand, M. Sigrist, and S. Haas, Phys. 
Rev. Lett. 86, 1086 (2001). 



